The magnetic fields of accreting T Tauri stars 
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Abstract. Models of magnetospheric accretion on to classical T Tauri stars often assume that the 
stellar magnetic field is a simple dipole. Recent Zeeman-Doppler imaging studies of V2129 Oph 
and BP Tau have shown however that their magnetic fields are more complex. V2129 Oph is a high 
mass T Tauri star and despite its young age is believed to have already developed a radiative core. In 
contrast to this, the lower mass BP Tau is likely to be completely convective. As the internal structure 
and therefore the magnetic field generation process is different in both stars, it is of particular interest 
to compare the structure of their magnetic fields obtained by field extrapolation from magnetic 
surface maps. We compare both field structures to mulitpole magnetic fields, and calculate the disk 
truncation radius for both systems. We find that by considering magnetic fields with a reaUstic 
degree of complexity, the disk is truncated at, or within, the radius obtained for dipole fields. 
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INTRODUCTION - MAPPING STELLAR MAGNETIC FIELDS 

Classical T Tauri stars (cTTs) are young pre-main sequence stars which accrete material 
from circumstellar disks. Many observations support the magnetospheric accretion sce- 
nario, where the stellar magnetic field is strong enough, and sufficiently well-ordered, 
to truncate the disk at a few stellar radii [IJ. Gas is forced to follow the field lines of the 
stellar magnetosphere and accretes on to the star at high velocity, producing detectable 
hotspots [2] . However, until recently little was known about the geometry of their mag- 
netic fields, with the majority of models assuming a dipole. 

Measuring the Zeeman broadening of unpolarized spectral lines has proved to be very 
successful in detecting T Tauri magnetic fields (e.g. [3]), although this technique is 
insensitive to the magnetic topology. In contrast, measuring the polarization signature 
in spectral lines gives access to the field topology, and allows information about how the 
magnetic energy is distributed within the different field components to be determined. 
However, like all polarization techniques, this suffers from flux cancellation effects and 
yields limited information regarding the field strength. 

The polarization signatures detected in spectral lines are typically small, and therefore 
cross-correlation techniques (such as Least-Squares Deconvolution [4J) were developed 
in order to extract information from as many spectral lines as possible. The signal-to- 
noise ratio of the resulting average Zeeman signature is several tens of times larger 
than that of a single spectral line. Magnetic surface features produce distortions in the 
Zeeman signature that depend on the latitude and longitude of the magnetic region, 
as well as the orientation of the field lines. By monitoring how such distortions move 
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FIGURE 1. Flattened polar projections of V2129 Oph (left) and BP Tau (right) showing the strength 
in Gauss and polarity of the radial field component |5 6|. Red denotes positive, and blue negative, field 
regions. The solid black circle represents the stellar equator, and the dashed circles lines of constant 
latitude. The surface magnetic fields are clearly non-dipolar. 



through the Zeeman signature as the star rotates, a method referred to as Zeeman- 
Doppler imaging (ZDI), the distribution of magnetic polarities across surface of cool 
stars can be determined (see Donati et al, these proceedings). 

ZDI studies have now been carried out on two cTTs, V2129 Oph and BP Tau [|5l|6l. 
Both stars have complex non-dipolar magnetic fields (Fig. [TJ. The field of V2129 Oph 
is dominated by a 1.2 kG octupole field component, tilted at ~ 20° with respect to the 
stellar rotation axis. The dipole component was found to be weak, with a polar strength 
of only 0.35 kG and tilted at ~ 30° IS. The field of BP Tau consists of strong dipole and 
octupole field components of strength 1.2 and 1.6 kG, both tilted by ~ 10° with respect 
to the stellar rotation axis, but in different planes [!6|. 



Field extrapolation and comparison with multipoles 

From the surface maps of the photospheric field, we reconstruct the coronal field 
via field extrapolation, using the potential field source surface method [7J. This has the 
advantages over full MHD models of computational speed and simplicity, and does not 
require assumptions about an equation of state, or the energetics. The disadvantages are 
that non-potential and time dependent effects can not be modeled. However, it has been 
demonstrated that the potential field model is often adequate and produces results which 
match those of more complex MHD models [O. We note that for cTTs the main source 
of non-potentiality in the field is likely to be due to the interaction between the stellar 
magnetosphere and the disk. We do not consider this here - the field structures that we 
present are static and do not evolve in time. 

The potential field source surface model requires two boundary conditions. The first 
is that the radial component of the magnetic field at the stellar surface is that which is 
determined from the Zeeman-Doppler maps. The second is that at some height above the 
star, known as the source surface Rs, the field becomes purely radial, which mimics the 
effect of a stellar wind blowing open field lines at some height above the star. Although 



FIGURE 2. The field structure of V2129 Oph (left) and BP Tau (right) derived by field extrapolation 
from the magnetic surface maps in Fig.[T| Closed/open field lines are shown in white/blue. There is a clear 
distinction between the complex surface field and the simpler larger scale field. 



the choice of Rs is a free parameter of the field extrapolation model, its location is well 
constrained by the observed accretion hotspot locations (Si O. Once Rs has been set, 
predictions can then be made regarding the X-ray emission properties of the coronal 
field [|9l, which will be tested against future X-ray observations. 

The field extrapolation technique applied specifically to cTTs is discussed in detail by 
[fTOl . while Fig. [2] shows the field topology obtained by extrapolation from the magnetic 
maps of V2129 Oph and BP Tau shown in Fig. [T| For both stars the surface field is 
complex and non-dipolar, while the larger scale field appears more well-ordered (Fig. 
[2]). As magneto spheric accretion models often assume that cTTs have simple dipolar 
fields (with the exception of recent work, e.g. [1 IJ), it is of interest to compare the field 
structures obtained by field extrapolation to a dipole with a source surface. Furthermore, 
as both stars have strong octupole field components, we also compare the extrapolated 
fields with octupoles. To compare the field structures we calculate the heights and widths 
of the field line loops obtained by field extrapolation (Fig. [3]). We define the height, h, 
of a field line loop as the maximum height of the loop above the stellar surface. The 
width of a field line, w, is defined as the distance along the segment of the great circle 
connecting the field line footpoints on the stellar surface [IT2II . 

For BP Tau the field structure follows a similar trend to a dipole field, which reflects 
its very strong dipole field component. However, the field structure of V2129 Oph shows 
significant departures from a dipole, and is characterized by numerous compact loops. 
The difference in field complexity is likely related to the different internal structure 
of both stars, which affects the magnetic field generation process, and ultimately the 
topology of the large scale field. V2129 Oph is believed to have a radiative core [5 J and 
has a particularly complex field structure, whereas BP Tau is completely convective and 
is found to have a simpler, although still non-dipolar, field topology [6J . 

For both V2129 Oph and BP Tau we find that larger scale field lines are wider than 
dipole magnetic loops (Fig. |3]). The larger scale loops, which are the ones that will 
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FIGURE 3. Height vs width of the field fine loops derived by extrapolation form the surface maps of 
V2129 Oph (left) and BP Tau (right). The black points represent the extrapolated fields and the red/green 
points a dipole/octupole fields. For both stars, the larger scale field lines are wider than a dipole. 



interact with the disk, are distorted close to the star by the strong surface field regions. 
Therefore, even though the large scale magnetic field is well-ordered and dipole-like, 
the structure of the accreting field is influenced by the complex field regions close to the 
star [[T2]| . We note that recently [13J have come (independently) to the same conclusion 
for some stars in their sample, and comment that in order to explain the deep absorption 
component detected in a certain line of He, "flows near the star with less curvature than 
a dipole trajectory seem to be required". 



Disk truncation 

In this section we compare how the disk truncation radius Rf differs between consider- 
ing a dipole magnetic field and the fields derived by field extrapolation. Various methods 
of calculating Rt have been developed [14], and in this work we assume that the inner 
disk is truncated where the torque due to viscous processes in the disk is comparable to 
the magnetic torque due to the stellar magnetosphere [fTSll . We assume that at the inner 
disk the perturbed toroidal component is equal to the poloidal component (B^ ~ B^, see 
[ll2|). This assumption should be valid as long as the disk is truncated sufficiently within 
the corotation radius and assuming the disk is strongly coupled to the stellar magnetic 
field [[T6J. Equating the differential magnetic and viscous torques, and assuming that the 
poloidal field threading the disk is dominated by the vertical component (B^ ^ B^), gives 

r^Bl = -M (^— ^ J . (1) 

The sohd lines in Fig. |4] represent the variation in the RHS equation ([T]) along the 
stars equatorial plane. We use the same stellar parameters as in [|T2ll . namely for 
BP Tau a mass, radius, rotation period and accretion rate of O.VM©, 1.95R0,7.6d 




FIGURE 4. The variation along the equatorial plane of the differential magnetic/viscous torque [i.e. the 
LHS and RHS of equation ^] for V2129 Oph (left) and BP Tau (right). The solid Une represents the RHS 
of ([T]), with the LHS for an aligned dipole field (dashed line), a tilted dipole field with a source surface 
(dash-dot line) and the field extrapolation of V2129 Oph/BP Tau (dotted line). In both cases the disk is 
truncated within the corotation radius of 6.7Rh. (V2129 Oph) and 7.4 R* (BP Tau). 



and 2.88 x lO^^Moyr^i respectively, and for V2129 Oph, 1.35M0,2.4R0, 6.53d and 
1 X lO^^Moyr-i. The various dashed/dotted lines in Fig. |4] represent the LHS of equa- 
tion ([T]) assuming different forms for the stellar magnetic field. The value of r where the 
lines cross, is then the disk truncation radius Rt. In the figures we have considered both 
an analytic dipole field, and a dipole field with a source surface surface, in order to allow 
a better comparison with the extrapolated fields of V2129 Oph and BP Tau. We assume 
that the dipoles have the same polar strength as the measured dipole components of each 
star (0.35/1.2 kG for V2129 Oph/BP Tau). 

For both V2129 Oph and BP Tau we find that the disk is truncated within the 
corotation radius, justifying our assumption above (Fig. |4]). For the extrapolated field 
of BP Tau, Rt is only slightly within what is expected for a dipole magnetic field. This 
reflects the strong dipole component of BP Tau, which is the dominant field component 
at the inner edge of the disk. In contrast, Rt for the extrapolated field of V2129 Oph is 
almost a stellar radii within that calculated using a dipole magnetic field. This is due to 
intrinsic complexity of the magnetic field of V2129 Oph. The field strength in the disk 
midplane for this star is less, as the drop-off in field strength with height above the star 
is greater than with BP Tau, and therefore the disk is able to extend closer to the stellar 
surface than would be expected with a dipole magnetic field [12J. Our conclusion is that 
when considering magnetic fields with a realistic degree of complexity the disk will be 
truncated at, or within, the radius determined from dipole magnetic field models. 



CONCLUSIONS AND FUTURE WORK 



We tentatively find that the complexity of the magnetic field topology of accreting T 
Tauri stars depends on the internal structure of the star. V2192 Oph, which has already 



developed a radiative core despite its young age, has a more complex and dominantly 
octupolar field, compared to the much simpler field of the completely convective BP Tau, 
which has a strong dipole field component. However, this conclusion is based on only 
two stars and more data across a range of stellar masses will be required to confirm if 
this is a general feature of all cTTs, although a similar result has been found for low- 
mass main-sequence stars (see Morin et al, these proceedings). For V2129 Oph there is a 
more rapid drop-off in field strength with height above the stellar surface compared with 
BP Tau, which means that the inner disk is truncated much closer to the star compared 
with a dipole magnetic field model. For BP Tau, with its stronger dipole component, Rt 
is roughly the same as expected with a dipole field. Thus the interaction between the 
stellar field and the disk may be different for stars with different internal structures. 

The surface fields of both stars are particularly complex, especially for V2129 Oph, 
but for both stars the larger scale, likely accreting field, is more well-ordered. However, 
the structure of the larger scale field is distorted close to the stellar surface by the strong 
surface field regions. Thus many of the larger scale field lines are wider than dipole 
magnetic loops, a result which other authors have recently found [fT3l . 

In order to quantify the above effects and discover if they are general features of all 
accreting T Tauri stars, more spectropolarimetric data of stars across a range of masses, 
radii, rotation periods and accretion rates are required. Such new data will be obtained 
over the coming years as part of the MaPP (Magnetic Protostars & Planets) project, a 
large program with the ESPaDOnS spectropolarimeter at the Canada-France-Hawai'i 
telescope (PI: J.-F. Donati), which will provide 690 hours of observing between late 
2008 and 2012. Such new data will provide an unparalleled opportunity to investigate 
magnetism on solar-like pre-main sequence stars, and give new insight into the history 
of our Sun and Solar System at a time when the planets were begining to form. 
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